arXiv: 1509.06416vl [astro-ph.SR] 21 Sep 2015 


Astronomy & Astrophysics manuscript no. CR'protostars'letter'arXiv © ESO 2015 

September 23, 2015 


Cosmic-ray acceleration in young protostars 

M. Padovani'’^, P. Hennebelle^, A. Marcowith^ & K. Ferriere"^ 

* Laboratoire Univers et Particules de Montpellier, UMR 5299 du CNRS, Universite de Montpellier, place E. Bataillon, cc072, 34095 
Montpellier, France 

e-mail: [Marco. Padovani , Alexandre. Marcowith] ©umontpellier. fr 
^ INAF-Osservatorio Astrofisico di Arcetri, Fargo E. Fermi 5, 50125 Firenze, Italy 
3 CEA, IRFU, SAp, Centre de Saclay, 91191 Gif-Sur-Yvette, France 
e-mail: heimebelle@cea.fr 

IRAP, Universite de Toulouse, CNRS, 9 avenue du Colonel Roche, BP 44346, F-31028 Toulouse Cedex 4, France 
e-mail: katia. ferriereSirap. omp. eu 

Preprint online version: September 23, 2015 


ABSTRACT 

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards 
the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent obser¬ 
vations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained 
by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be 
satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks 
can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is 
located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate 
protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can 
be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong 
and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules. 
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1. Introduction 

It is largely accepted that Galactic cosmic rays, which pervade 
the interstellar medium, are likely produced in shock waves 
in supernova remnants (Drury I1983I) . Cosmic rays activate the 
rich chemistry that is observed in a molecular cloud (Duley & 
Williams I19841 I and they also regulate its collapse timescale 
(Balbus & Hawley [1991 1 Padovani et al. l2014ll . determining the 
efficiency of star and planet formation, but they cannot penetrate 
in to the densest parts of a molecular cloud, where the formation 
of stars is expected, because of energy losses and magnetic field 
deflections (Padovani et al. 120091 Padovani & Galli l20lTl 120131 
Padovani et al. 120131 Cleeves et al. I2013I I. Recently, observa¬ 
tions towards young protostellar systems showed a surprisingly 
high value of the ionisation rate (Ceccarelli et al. 120141 Podio et 
al. 12014b . the main indicator of the presence of cosmic rays in 
molecular clouds. Synchrotron emission, the typical feature of 
relativistic electrons, was also detected towards the bow shock of 
a T Tauri star (Ainsworth et al. 12014b . Nevertheless, the origin of 
these signatures peculiar to energetic particles is still puzzling. 
Here we show that particle acceleration can be driven by shock 
waves occurring within protostars. 


increase, reaching the relativistic domain. The motion of par¬ 
ticles back and forth from upstream to downstream requires the 
presence of magnetic fluctuations that produce a scattering of 
the pitch angle, namely the angle between the particle’s velocity 
and the mean magnetic field (Drury [19831 Kirk ll994l l. We ar¬ 
gue below that the accelerated particles themselves can produce 
the necessary fluctuations to maintain DSA around shocks, as 
discussed, e.g., in Bell (11978b . 

In the following subsections we describe all the conditions 
that have to be satisfied in order to effectively accelerate pro¬ 
tons and electrons through DSA. The acceleration of helium and 
heavier nuclei will be presented in a following paper. All the 
constraints limiting the maximum energy of the accelerated par¬ 
ticles are written as functions of the upstream flow velocity in 
the shock reference frame, t/sh = - fsh, and t^sh are the 

flow and the shock velocities in the observer reference frame, 
respectively. Our aim is to compute the maximum energy, Emax, 
reached by a particle in the regime where ions and neutrals are 
coupled. As explained in Sect. 2.3 in this case the damping of 
the Alfven waves, which determines the confinement of parti¬ 
cles, is weak and particle acceleration is more effective. 


2. Particle acceleration in shocks 

This works focuses on shock acceleration by means of the diffu¬ 
sive shock acceleration (DSA) mechanism. Also known as first- 
order Fermi acceleration, DSA is a process where charged par¬ 
ticles systematically gain energy while crossing a shock front. 
Multiple shock crossings allow the particle energy to rapidly 


2.1. Condition on shock veiocity 

In order to have efficient particle acceleration, the flow has to be 
supersonic and super-Alfvenic. These two conditions are com¬ 
bined into the following relation 

t/sh,2 > max {9 X IQ-^ [y^l + , 2 x , 

( 1 ) 
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where t/sh ,2 has units of 10 ^ km s“', yad is the adiabatic index, 
Tn the upstream temperature in 10"^ K, ne the total number den¬ 
sity of hydrogen in 10 ^ cm“^, x the ionisation fraction, and 
the magnetic field strength in 10“^ G. The two terms in square 
brackets on the right-hand side of Eq. 0 are the ambient (or 
upstream) sound speed and the Alfven speed of the total gas in 
10 ^ km s“', respectively. 


2.2. Condition on low-energy particle acceleration: collisional 
losses 


We are interested in the acceleration of low-energy particles (< 
100 MeV-1 GeV), because they are responsible for the bulk of 
the ionisation. We have to verify that the shock acceleration rate 
is larger than the collisional loss rate Following 

Drury et al. (11996b . the acceleration rate is given by 


3.2x10-* k“(r-l) 

y-1 r(lH-rkd/ku) 




sh,2 


B 5 s- 


( 2 ) 


where fi - mlnip is the particle mass normalised to the proton 
mass, and are the diffusion coefficients in the upstream and 
downstream media, respectively, normalised to the Bohm value 
for protons 
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with y the Lorentz factor, p - y'* - 1, and r is the shock 

compression ratio. For a parallel shock, a - and 
while for a perpendicular shock, a = 1 and /Cu = r/Cd[^The energy 
loss rate is given by 


Cs = 3.2 X 10- 


>0 


fl n6T-25 S 


(4) 


where L _25 is the energy loss function (Padovani et al. 12009b 
in units of 10^^^ GeV cm^ which was extended to lower en¬ 
ergies including Coulomb losses (for protons, Mannheim & 
Schlickeiser 1 1 994| and electrons, Swartz et al. 11971b and syn¬ 
chrotron losses (Schlickeiser 12002b . The maximum energy of 
accelerated particles set by energy losses, Fioss, is found when 
face = specifically 


/3T-25 = 10 


K(r-l) 

r(l + rkdlkyP 




(5) 


2.3. Condition on particle acceleration: ion-neutral friction 

The main limit on the possibility of particle acceleration is given 
by the presence of an incomplete ionised medium. In fact, the 
collision rate between ions and neutrals can be as high as to 
decrease the effectiveness of the DSA, damping the particle’s 
self-generated Alfven waves responsible for the particle scatter¬ 
ing back and forth the shock. Ions and neutrals are effectively 
decoupled if the wave frequency is larger than the ion-neutral 
collision frequency. Following Drury et al. (11996b and account¬ 
ing for the fact that particles are not fully relativistic, we find that 
the critical energy separating these two regimes, Fcoup^ is derived 
by solving the following relation 

yP = 8.5 X . ( 6 ) 

* A parallel/perpendicular shock is when the shock normal is paral¬ 
lel/perpendicular, respectively, to the ambient magnetic field. 


If the particle energy is larger than Fcoup, ions and neutrals are 
coupled. 

The upper cut-off energy due to wave damping, Fdamp, is set 
by requiring that the flux of accelerated particles advected down¬ 
stream by the flow is equal to the flux of particles lost upstream 
because of the lack of waves (due to wave damping) to confine 
the particles. Following Drury et al. (11996b . using their exact 
equation for the wave damping rate, accounting for departures 
from fully relativistic behaviour, and assuming U^h to be much 
larger than the Alfven speed, Fdamp follows from 

yP^ = 8.8 X - x)-^B-_\P'_^ , (7) 

where 

S = BI 5 -H 1.4 X IQ^^p^y^p^T^-^nl?? , ( 8 ) 

Both Eq. ^ and Eq. § are valid for T e [10^, 10^] K. P'_^ 
is the fraction of the shock energy {mptiU'^^ going into particle 
acceleration in 10^^. P' is proportional to the shock efficiency 
77 e [10 ®, 10 *] (Bykov l2004b . which represents the fraction of 
particles extracted from the thermal plasma and injected into the 
acceleration process by a shock. We predict both non-relativistic 
and mildly relativistic accelerated particles and we checked a 
posteriori that there is no strong back-reaction. This means that 
the upstream medium is not warned by these particles that a 
shock is coming and we can safely assume that the shock and 
the DSA process are unmodified. In other words, calculations 
are carried out in the test-particle limit. 

If Fdamp > Fcoup, then Fdamp IS in the coupled regime, namely 
neutrals coherently move with ions and ion-generated waves are 
weakly damped. The last inequality can be written by combining 
Eqs. (j^ and (0 as 

10 * 

^ = —Ht/^\ 2 « 6 X' ®(1 - x)-'b:®p :2 > 1. (9) 

We consider shocks in three types of environments; in jets as 
well as in accretion flows in the collapsing envelopes and on the 
surfaces of protostars. Using the range of parameters of Table[T] 
we find that 1 in protostellar envelopes (Sect.|^. This is to 
say that the following two conditions on shock age and geometry 
(Sect. [ 2 ^ are only discussed with reference to shocks in jets and 
on protostellar surfaces. 


2.4. Conditions due to shock age and geometry 

The maximum energy set by the age of the shock, Fage, is found 
when the acceleration time, given by the inverse of Eq. (|^, is 
equal to the age of the shock. The latter can be assumed of the 
order of the dynamical time of the jet (> 10* yr, de Gouveia 
Dal Pino 11995b or equal to the accretion time in the case of a 
surface shock (~ 10* yr, Masunaga & Inutsuka 12000b . Then, 
Fage is computed from 


2 k“(r-l) , 2 


( 10 ) 


with fage 3 in units of 10 * yr. 

A further constraint is given by the geometry of the shock. 
In particular, the upstream diffusion length, Tu = tCu/Ush, has 
to be at most a given fraction e < 1 of the shock radius, Rsh', 
besides, in the jet configuration particles may also escape in the 
transverse direction (the shock can be assumed planar as long 
as the particle’s mean free path around the shock is smaller than 
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the transverse size of the jet, Rj_). The maximum energy due to 
upstream escape losses, fiesc.u, follows from 

yj3^ = 4M“il^^Ush.2B-5 min (e/?sh, 2 ,^±, 2 ) , ( 11 ) 

where both Rsh ,2 Rj^ 2 in units of 10^ AU. In the follow¬ 
ing we assume e = 0.1 (Berezhko et al. 11996b . Since jet shocks 
have a small transverse dimension, there is a further condition 
for the escape of particles downstream: the maximum energy 
due to downstream escape losses, itesc.d, is found when the ac¬ 
celeration time, inverse of Eq. @, is equal to the downstream 
diffusion time, ^ - /?^/(4A'd)Qnamely 

y/3\y - 1) = 5.8 ~^h -HU,h, 2 B- 5 R^. 2 f ■ (12) 

r(l H- rkd/ku) 

Finally, if the shock is supersonic and super-Alfvenic (Eq.[T]l and 
if^ > 1 (Eq.|g), the maximum energy reached by a particle is 
Emax “ min[£'ioss, Edajnp? Eage, Egsc.u? Egsc.d]- 

3. Potential particle acceleration sites 

In this Section, we identify and characterise possible sites of par¬ 
ticle acceleration in protostars. In particular, we consider accre¬ 
tion flows in the collapsing envelopes and on protostellar sur¬ 
faces as well as jets. The required parameters needed to prove 
the effectiveness of shock acceleration are shown in Table [T] 

It is straightforward to verify that Eq. (|^ is not fulhlled in en¬ 
velopes «c 1). The ionisation fraction and the shock velocity 
are too small, quenching the particle acceleration. The magnetic 
held strength is also as large as to produce a sub-Alfvenic shock. 
This rules out envelopes as a possible shock acceleration sites. 
From now on we will focus on shocks in jets and on protostellar 
surfaces. 


Table 1. Values of the parameters described in the text. 


site* 


T 

n 


X 


B 


[km s“*] 

[K] 

[cm" 




[G] 

£ 

1 - 10 

50-100 

10' - 

1(F 

< 10 

-6 

1 

0 

1 

T 

0 

J 

40-160 

10“* - 10^ 

10" - 

10^ 

0.01 - 

0.9 

5 X 10-5 _ ^0-3 

S 

260 

9.4 X 10^ 

1.9 X 

10‘2 

0.01 - 

0.9 

0 

1 

1 

0 


*£ = envelope; = jet; <S = protostellar surface. 


3.1. Jets 

Jets are observed at all stages during the evolution of a protostar 
(e.g. McCaughrean et al. I20021 Reipurth et al. ll997l Watson & 
Star)elfeldt l2004b . Jet speeds, tJfl, are similar for different classes 
(60 - 300 km s"') with shock velocities, Ush, of the order of 
20-140 km s ' (Raga et al.l2002l 120141 Hartigan & Morse l2007l 
Agra-Amboage et al. 1201 lb . then t/sh = 40 - 160 km s '. The 
neutral density is between 10"' and 10^ cm~^ (Lehoch et al. l20121 
Gomez-Ruiz et al. 12012b with temperatures of the order of 10"' K 
up to about 10® K (Frank et al. 12014b . There is only one theoret¬ 
ical estimate for the magnetic held strength (B ~ 300 - 500 fiG) 
for Class II sources (Te§ileanu et al. 12009112012b . The transverse 
radius of a jet is about 5 AU, 10 AU, and 50 AU at 15 AU, 
100 AU, and 1000 AU from the source, respectively (Cabrit et 
al. 120071 Hartigan et al. 12004b . The ionisation fraction in Class I 
and II are similar, x ~ 0.05 - 0.9 (Nisini et al. 120051 Maurri et 
al. 12014b . while Class 0 jets are mainly molecular (Dionatos et 
al. l20T0b . 

^ The factor 4 in the denominator comes from the fact that the diffu¬ 
sion process in the perpendicular direction is in two dimensions. 


3.2. Accretion flows on protostellar surfaces 

We use the computational results of the protostellar collapse of 
an initially homogeneous cloud core described in Masunaga & 
Inutsuka (12000b . Their simulation describes the phase of main 
accretion, when the protostar mass grows because of the steady 
accretion from the infalling envelope. They give the temporal 
evolution of temperature, density, and How velocity which, as¬ 
suming a stationary shock, is equal to the shock velocity. The 
radius of the protostar is set to 2 x 10“^ AU and we hnd that only 
the last time step of the simulation, corresponding to the end of 
the main accretion phase, leads to a strong proton acceleration. 
Parameters are listed in the third row of Table [T] 


4. Maximum energy of the accelerated particles 

For jet shocks we perform a parametric study using the val¬ 
ues in the second row of Table assuming a parallel shock, 
77 = 10 and T — 10"' K. We also consider - tcb, which 
is the most favourable circumstance for accelerating particles 
(Drury et al. 11983b and we compute E^ax for R^h = 100 AU 
and Rx = 10 AU. The upper panel of Fig. [T] shows the maxi¬ 
mum energy which a shock-accelerated proton can reach. By in¬ 
creasing both t/sh and B, E^ax attains higher values up to about 
13 GeV for protons. Once the combination of parameters satis- 
hes the condition > 1, E^ax rapidly reaches a constant value, 
encompassed by the cyan contours in each subplot. In fact, the 
maximum energy is controlled by Eesc.d that is independent of 
both n and x. Supposing the magnetic held to have a strong 
toroidal component, we repeat the calculation for the case of a 
perpendicular shock hnding that E^ax decreases by a factor of 
about 1.3, since oc (r + l)/(2r)j^ 

For shocks on protostellar surfaces, we use values in the 
third row of Table [l] varying x and B. Assuming 77 = 10^®, 
ku = 1 , and a parallel shock, we hnd values of Emax for pro¬ 
tons up to about 26 GeV for Z? ~ 3 - 10 G (see lower panel 
of Fig. [T]i, which are comparable with magnetic held intensities 
computed by e.g. Garcia et al. (12001b . Because of high tempera¬ 
tures, Coulomb losses are dominant and E^ax is constrained by 
Eioss. Thus, for a perpendicular shock, E^ax is a factor of about 
1.3 larger, since cx (r -(- l)/r. Hatched areas in both 

panels of Fig. show regions where acceleration is not possible 
because of strong wave damping. 

Electrons can be accelerated as well, but generally E^ax for 
electrons is much smaller than E^ax for protons because of wave 
damping and stronger energy losses. As an instance, for t/sh = 
160 km s ' and B - 1 mG, Emax ~ 300 MeV for a narrow range 
of density and ionisation fraction (77 > 3 x 10® cm“^, x > 0.6). 
For lower values of B and t/sh, Emax ^ 50 MeV. We also hnd 
that electron acceleration is not triggered by protostellar surface 
shocks. 


5. Discussion and conclusions 

We investigated the possibility of accelerating particles within 
a protostellar source by means of shock processes through the 
diffusive shock acceleration mechanism. We focused our atten¬ 
tion on the effectiveness of shocks in envelopes, on protostellar 
surfaces, and in jets. We concluded that: 


^ Superscripts _L and || refer to perpendicular and parallel shocks, re¬ 
spectively. 
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Fig. 1. Upper panel. Case of a parallel shock in jets: ioni¬ 
sation fraction, x, versus density of neutrals, n, for t/jh = 
40 and 160 km s ', B = 50 pG and 1 mG, and T - 10"^ K. Lower 
panel. Case of parallel (left) and perpendicular (right) shocks 
on protostellar surfaces: ionisation fraction, x, versus magnetic 
field strength, B, for U^h = 260 km s ', T = 9.4 x 10^ K, 
and n - 1.9 x 10'^ cm“^. Colour maps show values of Fmax 
for protons and cyan contours delimit the regions where Fmax 
reaches its maximum constant value shown in GeV in each sub¬ 
plot. Vertically hatched regions refer to combinations of param¬ 
eters corresponding to strong wave damping (1% < 1). 


(i) In envelopes, x and Ush are too small, preventing any particle 
acceleration. Besides B is as large as to yield sub-Alfvenic 
shocks. 

(ii) Jet shocks are possible accelerators of particles that can be 
easily boosted up to relativistic energies. The acceleration is 
more efficient for protons which can reach up to about 13 or 
10 GeV (for parallel or perpendicular shocks, respectively), 
while electrons attain at most about 300 MeV because of 
wave damping and energy losses. 

(Hi) Protostellar surface shocks can accelerate protons up to 
about 26 or 37 GeV (for parallel or perpendicular shocks, re¬ 
spectively). Electrons cannot be accelerated mainly because 
of large magnetic field strengths leading to synchrotron 
losses. 


The set of conditions that has to be fulfilled is highly non¬ 
linear: small variations in one or more parameters (B, x, n, T, 
Ush, n, ki) can make the acceleration process inefficient. As 
a consequence, since young protostars are highly dynamic sys¬ 
tems, particle acceleration can be a very intermittent process. In 
a following paper we will discuss in detail other possible ac¬ 
celeration mechanisms as well as the effect of variations in the 
parameter set and the distance of the jet shocks from the proto¬ 
star, including departures from Bohm-type diffusion with a view 
to studying the propagation of high-energy particles in the pro¬ 
tostellar environment in order to explain the available observa¬ 
tions. 
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